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Following erratum submitted 29th January: 



While reading the description of Figure 2, as well as the caption of the same 
Figure, it is to be noted that the solid, dashed, short dashed and dotted lines were 
drawn for = 0.001, 0.01, 0.1 and 1.0 respectively. The sequence in inadvertantly 
reversed in the text. 
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Abstract 

We study the spectral properties of a very general class of accretion disks 
which can be decomposed into three distinct components apart from a shock at 
r = r s : (1) An optically thick Keplerian disk on the equatorial plane (r > r s ), (2) 
A sub-Keplerian optically thin halo above and below this Keplerian disk r > r s 
and (3) A hot, optically slim, r ~ 1 postshock region r < r s ~ 5 — Wr g where r g 
is the Schwarzschild radius. The postshock region intercepts soft photons from 
the Keplerian component and reradiates them as hard X-rays and 7 rays after 
Comptonization. We solve two-temperature equations in the postshock region 
with Coulomb energy exchange between protons and electrons, and incorporat- 
ing radiative processes such as bremsstrahlung and Comptonization. We also 
present the exact prescription to compute the reflection of the hard X-rays from 
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the cool disk. We produce radiated spectra from both the disk components as 
functions of the accretion rates and compare them with the spectra of galactic 
and extragalactic black hole candidates. We find that the transition from hard 
state to soft state is smoothly initiated by a single parameter, namely, the mass 
accretion rate of the disk. In the soft state, when the postshock region is very 
optically thick and cooled down, bulk motion of the converging flow determines 
the spectral index to be about 1.5 in agreement with observations. 

Subject headings: accretion, accretion disks - black hole physics - radiation mech- 
anisms: nonthermal - shock waves -stars:neutron 
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1. INTRODUCTION 



In the 1970s, the standard accretion disk models were constructed (Shakura & 
Sunyaev 1973; Novikov & Thorne 1973) to largely explain observations of binary sys- 
tems. Matter supply from the companion is assumed to be through the Roche lobe and 
the angular momentum distribution is assumed to be Keplerian throughout the disk. 
Here, viscosity drives the inflow by removing angular momentum outward and keeping 
the entire disk Keplerian in the process. In the early 1980s, these disks became favorite 
candidates for the explanation of the "big blue bump" seen in the UV region of the 
active galaxies (Malkan 1982; Sun & Malkan 1989 and references therein) with some 
success. 

With the advent of space based observations, the X-ray and 7-ray spectra of 
these objects show that there could be even bigger bumps in the X-ray regions of the 
spectra. Galactic black hole candidates such as GS1124-68, show two distinct spectral 
components, namely, soft and hard, which apparently vary quite independently (Tanaka 
1989, hereafter TAN; Ebisawa et al, 1993, 1994 hereafter E93 and E94, respectively; 
Inoue 1992, hereafter 192). It was becoming clear that the standard accretion disks, 
which are in general cooler, cannot explain the production of X-rays and 7-rays. These 
models also cannot explain the almost zero time-lag correlated variabilities observed in, 
e.g., NGC4151, NGC5548 (Clavel et al. 1990; Peterson et al. 1991; Perola et al. 1986). 
Explanations of the hard components seen in the spectra of neutron star candidates 
and galactic and extragalactic black hole candidates required extra components, such 
as plasma clouds, hot coronae etc. which are not self-consistently constructed. 

That the disks need not be of "standard" type was sensed by theoreticians since 
early 1980s (Liang & Thomson, 1980; Paczyhski & Bisnovatyi-Kogan 1981; Muchotrzeb 
& Paczyhski 1982; Paczyhski & Wiita 1980; Abramowicz et al. 1988, Chakrabarti 1989, 
hereafter C89; Chakrabarti 1990ab, hereafter C90ab, Chakrabarti & Molteni 1995, here- 
after CM95). Particularly, tt was realized that the disks are not Keplerian everywhere, 
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especially at the inner boundary. (For a general discussion see, Chakrabarti 1995a, 
hereafter C95a). First global and fully self-consistent solution of viscous, transonic 
flows was obtained by C90ab who found that a Keplerian flow at the outer bound- 
ary can become sub-Keplerian at the inner edge (with or without shock waves). Or, 
contrarily, an initially sub-Keplerian flow can become almost Keplerian when viscosity 
is high enough. In a binary system, which includes a black hole or a neutron star, 
matter could be accreted from the winds of the companion at the same time as it is 
through the Roche lobe. In low mass X-ray binaries (LMXB), where companion winds 
may be scarce (except probably during outburst phase) the sub-Keplerian flow could 
come from the Keplerian disk itself as the advection becomes important closer to the 
black hole. In an active galaxy, the incoming matter is either a mixture of Keplerian 
and sub-Keplerian flows, or, entirely sub-Keplerian. This is because matter may be 
supplied from constantly colliding winds from a large number of stars ~ pc away and 
thus most of the angular momentum of the infalling matter may be lost. The above 
solutions, obtained in the context of isothermal flows, have been generalized since then 
for any heating and cooling (C95a; Chakrabarti 1995b, hereafter C95b) but no new 
topologies emerged. 

After the matter starts with sub-Keplerian or Keplerian angular momentum I < 
lKep(r ut), its subsequent behavior depends strongly upon the accretion rate and the 
viscosity in the flow (C90ab; CM95, C95ab, Appendix). A unified description of what 
the flow might do (e.g., Fig. 4 of C89) suggests that, for a given entropy of the flow, 
if the accretion rate is small enough, the flow passes through the outer sonic point 
(just as in a Bondi flow) and remains supersonic before falling onto a black hole. If 
the accretion rate is high, matter would pass through the inner sonic point. In the 
intermediate case, matter is virtually stopped at a distance r S) where / — lKe P (r s ) due 
to the centrifugal barrier, and a shock is formed. The postshock flow becomes highly 
subsonic momentarily, but picks up its velocity as it approaches the black hole eventually 
entering the horizon supersonically. The postshock region behaves similar to the thick 
accretion disks discussed in the literature (e.g., Paczyhski & Wiita 1980) but is more self- 
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consistent as advection is explicitly added (C89, Molteni, Lanzafame, & Chakrabarti 
1994, hereafter MLC94). In these solutions, the entropy of the flow passing through the 
inner sonic point is higher compared with the entropy of the flow passing through the 
outer sonic point. The combined effects of the shock and the viscosity generate the right 
amount of entropy to allow the flow to pass through the inner sonic point. The energy 
of the flow is either constant or very nearly so, because the energy and the entropy are 
advected with the flow toward the black hole. The advection is clearly important, since 
the radial velocity could be very high close to the black hole. When accretion rate is 
low enough, and shock conditions are not satisfied, a Keplerian disk would become sub- 
Keplerian close to the black hole at the same time becoming gas pressure dominated 
(Rees et al., 1982). Recently, importance of advection is stressed ('newly discovered 
advection dominated flows') by a simple self-similar analysis by Narayan & Yi (1994) 
who analyzes one of the branches of the solution topologies found by C90a,b. 

The most general and complete global solution to-date of the viscous sub-Keplerian 
flow with advection (which becomes Keplerian at a larger distance), shows a very inter- 
esting behavior (C90a,b, CM95, C95b). When the viscosity is low [typically, when the 
viscosity parameter (Shakura & Sunyaev 1973) is less than a critical value, a v < a vc ], 
the shock can still form (depending on accretion rate, see C89) as in an inviscid flow, 
but as the viscosity in increased (a v > a vc ), the stable shock becomes weaker until it 
finally disappears. Numerical simulation provides insight into what actually happens 
(CM95). In the presence of viscosity, the rate of transport of angular momentum is 
faster in the postshock subsonic flow as compared with the preshock supersonic flow. 
As a result, angular momentum starts 'piling up' in the postshock flow and pushes the 
shock away to a larger distance. If one employs the Shakura- Sunyaev (1973) prescrip- 
tion, and the viscosity is small enough (a v £ a vc ~ 0.01 for isothermal disk, C90a), the 
transport rates could be matched and the shock settles down farther out. If the viscos- 
ity is higher, the transport rate is fast enough to push the shock to extinction and the 
flow becomes subsonic and Keplerian except close to the inner edge of the disk, where 
it passes through the inner sonic point. In the case of more general cooling and heating 
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laws, a vc could be anywhere between to 1 depending upon the relative importance of 
cooling and heating (C95a,b). The predictions of shock locations and other properties 
from our vertically averaged model were found to agree with three dimensional simu- 
lations (MLC94). We, therefore, believe that these solutions are generally applicable 
even for quasi-spherical flow. This has also been noticed by recent self-similar study of 
Narayan & Yi (1994). 

Numerical simulations (MLC94, see also Hawley et al. 1984, 1985; Ryu et al. 
1995) also show that a strong wind forms in the postshock flow. Analytical works 
(C89) show that winds form when the specific energy of matter is positive and the 
specific entropy is high. The wind is found to become supersonic within a finite distance 
(MLC94). In the unified description (Fig. 4 of C89), the region of parameter space 
which produce winds (with and without shock waves) is discussed as well. 

These, exact and global solutions (C90a,b, C95a,b) with higher viscosity, of the 
original model of transonic viscous disks of Paczyhski and collaborators represent the 
most complete, Keplerian and quasi-Keplerian accretion disks in the literature. It has 
been suggested (Chakrabarti, 1993, hereafter C93, Chakrabarti, 1994, hereafter C94, 
C95a,b) that they go over to the so-called "thin" , "thick" , "slim" and "newly discovered 
advection dominated" and "cooling dominated" disks when appropriate approximations 
are made. Moreover, the solution of transonic viscous disk (C90a,b, CM95, C95a,b) also 
includes stationary shock waves for lower viscosity and appropriate accretion rate limit 
(C89). This 'grand unification' of disk models was possible, because the most general 
equations with all the solution branches including shocks are used (C93, C94, C95a,b) 
without making any further restrictive assumptions such as self-similarity (Narayan & 
Yi, 1994) or Keplerian distribution (Chen et al., 1995). Based on these solutions, one 
can construct the most general accretion disk close to a black hole. Assume for the 
moment that the flow is sub-Keplerian at the outer boundary and assume that the 
viscosity is higher on the equatorial plane and gradually drops with vertical height. 
The higher viscosity in the equatorial region re-distributes the sub-Keplerian angular 
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momentum and makes a standard, optically thick, Keplerian disk on the equatorial 
plane (CM95). However, the smaller viscosity at higher latitudes would produce an 
optically thin halo and this gas will pass through a shock close to the black hole, 
forming a very hot corona. In cases when an admixture of Keplerian and sub-Keplerian 
matter is supplied at the outer boundary (e.g., when both Roche lobe accretion and 
wind accretion are present) the segregation of this flow into two components is trivially 
achievable. In case where only Keplerian matter is supplied at the outer boundary, the 
disk can still be segregated in the same manner closer to the black hole depending on 
the viscosity variation and its vertical distribution which determines the region from 
which the flow starts deviating from Keplerian (C90a,b, Appendix). The solution with 
the shock (if available) is always chosen in place of the shock-free solution, as the former 
has a higher entropy (C89, MLC94). 

Fig. 1 shows a cartoon diagram of the generalized accretion disk in presence of 
Keplerian and sub-Keplerian matter close to a black hole (also see, C94). An optically 
thick Keplerian disk is flanked by an optically thin, sub-Keplerian halo which passes 
through a standing shock close to the centrifugal barrier. The postshock region forms 
a hot ion pressure supported torus (e.g. Rees et al., 1982; MLC94) of moderate optical 
depth t ~ 1 as the radial kinetic energy of the infall is thermalized: pv^(r s ) ~ P ~ 
nnkTp. For matter with marginally bound angular momentum, the shock forms around 
r s ~ 10r g (r g = 2GM/c 2 , the radius of the black hole) in Schwarzschild geometry. 
The flow remains dissipative up to the marginally stable radius r ms = 3r g , beyond 
which the flow rapidly passes through the sonic point. For a flow around a Kerr black 
hole of angular momentum parameter a = 0.99, these length scales are roughly half as 
long. (Henceforth, we choose units of distance and velocity to be r g and c respectively. 
We continue to use r to denote radial distances with this unit.) For higher angular 
momenta, shocks form farther away (C89). The postshock flow becomes geometrically 
thick, intercepting a few ( ;$ 5) percent of the soft photons from the optically thick, 
Keplerian component of the preshock accretion disk. The postshock region with a 
Thomson optical depth of r h ~ OApdx (~ 1 — 3 for M ~ M Edd rates) heats up 
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these soft photons through inverse Compton scattering (Sunyaev & Titarchuk 1980, 
hereafter ST80; Titarchuk 1994, hereafter T94; Titarchuk & Lyubarskij 1995; hereafter 
TL95 respectively.) and the photons are re-emitted which we claim to be responsible 
for the hard radiations from the accretion disks around galactic and extragalactic black 
holes. Some of these hard photons are intercepted by the Keplerian disk, a part of 
which is absorbed, thermalized, and reradiated by the disk as soft radiations, while 
the rest is reflected. We compute, completely self-consistently, the fraction of incoming 
hard radiation that each annulus absorbs or reflects away. Based upon our analysis, 
we do not find the component of hard radiation reflected from the disk to contribute 
significantly to the observed spectra. Secondly, we find the Keplerian disk component 
is not suitable for the site of the iron line emissions. Rather, we believe that the 
outflowing winds from the postshock flow (MLC94) or the evaporated preshock disk, 
could be the site of these lines at various stages of ionization (Chakrabarti et al., 
1995, hereafter CTKE95). We shall touch upon this aspect of the problem in the final 
section. When the net accretion rate is high enough, cooler photons trapped in our 
postshock region are Comptonized by the 'converging flow' (e.g. Blandford & Payne, 
1981) for r < r convg ~ 3(M d + M h )/M Edd r g where, M d and M h are the disk and halo 
rates respectively, and Mem denotes the Eddington accretion rate. The resulting hard 
spectrum is very weak with a spectral slope of a ~ 1.5 [F v oc v~ a ) and is characteristics 
of the spherical converging optically thick flow (Titarchuk, Mastichiadis & Kylafis 1995, 
hereafter TMK96). 

The advantage of this model is that the soft and the hard radiations are formed 
self-consistently from the same accretion disk without invoking any ad hoc components 
such as the plasma cloud, hot corona, etc., whose origins have never been made clear. 
Secondly, our results are sensitive to only one parameter, namely, the disk accretion rate 
m d , although some observations suggest a variation of the halo accretion rate m h cannot 
be ruled out. In the next Section, we present the equations for the hydrodynamics and 
radiative processes which we use to describe our model. In §3, we present the results 
and physical interpretations of our solution. In §4, we present general comparisons with 
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observations. In §5, we make concluding remarks. 

2. THE MODEL DESCRIPTION AND THE BASIC EQUA- 
TIONS 

2.1 Description of the Model 

As mentioned in the introduction, our model consists of two major disk compo- 
nents: 

(1) The Standard Disk Component: A standard, optically thick disk could be produced 
either from the Keplerian or the sub-Keplerian matter at the outer boundary, if the 
viscosity and cooling are appropriate (C90a,b, CM95, C95a,b, Appendix). While con- 
structing the equatorial Keplerian disk, we therefore incorporate the usual assumptions 
regarding the flux emitted from the disk. We assume that before the shock is formed 
at r = r s , there is a thin, standard, Keplerian disk on the equatorial plane which, in the 
absence of the shock, would emit a flux of (e.g, Shapiro & Teukolsky 1984) 

F SS = 7.6 x u^^^f—M—^ cm- s"'. (1) 

Here 1=1- (3/r) 1 / 2 . In the presence of the shock at r s , this component continues 
till r s . The soft radiation luminosity Ls is then the integral of the above flux from r s 
outwards plus the intercepted radiation coming from the postshock region: 

L s = L ss + L H f sd (l-A u ), (2) 

where, A v is the albedo of the disk, L H ~ L S sfdsEcomp is the hard luminosity from 
the postshock region, f s d is the fraction (determined by the geometry) of Lh that is 
intercepted by the disk, and Ecomp is the factor by which the soft photons incident on the 
postshock region (fraction f ds of L S s) is enhanced due to Comptonization. The second 
term therefore represents the luminosity of the soft photons resulting from reprocessing 
of the hard radiation intercepted by the Keplerian disk. 
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(2) Sub-Keplerian Halo Component: We consider in detail the properties of the sub- 
Keplerian halo component and the effects of the halo on the emission properties of the 
Keplerian component. We assume the halo to be quasi-spherical and axisymmetric. For 
the purpose of obtaining the height of the shock, we assumed the flow to be in vertical 
equilibrium in the postshock region. The height of the shock enables us to compute 
the fraction f S( i of disk soft photons intercepted by the shock as well as the fraction fd s 
of hard radiations from the shock intercepted by the disk. We assume the polytropic 
index to be 7 = 5/3 in the halo. A 7 = 5/3 spherical adiabatic flow cannot have a shock 
(C90b), but if the flow is thin, it can have two sonic points and therefore a standing 
shock. Even when a flow has only the inner sonic point, a shock can form if it is already 
supersonic at the source. An effective polytropic index is likely to be 4/3 > 7 > 5/3, 
but we do not expect any change of our general conclusions. Since we are interested in 
accreting solutions, the angular momentum in the halo component should be around the 
marginally bound and marginally stable values. In this low angular momentum regime, 
in the Schwarzschild geometry, the shock may be located around r s = 10—30 or so (C89). 
The inner edge of the postshock flow is extended till only the inner sonic point r iy since 
for r < Tj, advection is so strong that matter may be considered radiatively inefficient. 
Since the inner sonic point is close to the marginally stable orbit r ms = 3 anyway, we 
integrate the equations only from r s to r rns . We also considered the solutions around a 
Kerr black hole, where r s ~ 5 and r ms ~ 1.5 will be employed. These are reasonable 
values for a ~ 0.99. The results are similar, though vary only in details. Flow within 
the inner sonic point constitutes the so-called convergent inflow provided the accretion 
rate is high enough which produces the characteristic weaker hard component as will 
be discussed later. 



2.2 Hydrodynamics 



We solve two temperature equations for the sub-Keplerian halo component and 
the postshock region. We ignore the vertical structure of the halo and use the verti- 
cally integrated equations. Comparison of analytical (C89) and numerical simulations 
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(MLC94) show that such assumptions are justified. The radial momentum equation is 
given by, 



where P and p are vertically integrated pressure and density, respectively, / is the 
specific angular momentum and <fr(r) = — l/2(r — l)" 1 denotes the Paczyhski-Wiita 
(Paczyhsi & Wiita, 1980) potential which mimics the geometry around a black hole 
quite well. The variation of the azimuthal angular momentum in presence of viscosity 
is given by 



We assume the viscosity in the halo to be small enough (a v < a vc ) so that the centrifu- 
gal barrier can form a shock in the halo (C90a,b, CM95, C95a,b). However, we assume 
the viscosity on the equatorial plane to be high enough so that the angular momentum 
there is Keplerian. This would produce a generic composite disk with optically thick 
(Keplerian disk), thin (sub-Keplerian halo) and slim (r ~ 1, postshock region) compo- 
nents. Whether all the three components will be present will depend on viscosity and 
cooling processes. Detailed classification of global solutions will be discussed elsewhere 



In the region, r < r s , the Keplerian and sub-Keplerian components can behave 
identically for three reasons: (1) the angular momenta of the sub-Keplerian and Kep- 
lerian components are roughly the same at the centrifugal barrier, i.e., atr^r s ; (2) 
v r ~ in the immediately postshock region of the halo whereas dv r /dr in the 
disk component (velocity shoots up from a minimum, C89, C95b, Case A in Fig. 10 of 
Appendix A) — in either case the advection term v r dv r /dr, in the immediate vicinity 
of the shock is negligible and both the flows accrete downstream as a single component, 
and, most importantly, (3) since the postshock region is much hotter (more than 2 to 
3 orders of magnitude, see, §3) than the Keplerian component, it evaporates the disk 
underneath and together they are expected to behave as a single component (Fig. 1) 
provided the mixing of the disk and halo are efficient. In the present paper we assume 




(3a) 




(36) 



(C95b). 
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that the turbulent mixing time is slower than the infall time scale. In a neutron star 
accretion the mixing is complete. This case will be discussed elsewhere. 



The optically thin preshock halo does not radiate efficiently and therefore advects 
the energy and entropy along with the flow exactly as in our previous model of transonic 
disks (C89; C90a,b; C95b). In the postshock region, the disk will be hotter. If the flow 
is in hydrostatic equilibrium in the vertical direction, Mach numbers at the shock are 
related by (C89, C90b): 

[M + (3 7 -l) + 1 ^] 2 = [M_(3 7 -l) + 1 ^] 2 

2 + (7 - 1)M 2 2 + (7 - l)M 2 _ ' { a) 

(The subscripts "-" and "+" denote the preshock and postshock quantities, respec- 
tively.) Assuming a strong shock, it is easy to show that the post shock Mach number 
is M + ~ 0.35 which gives the postshock proton temperature as: 

T p ~ 1.56 x lO 12 ^? 1 K. (46) 

However, if the flow were not assumed to be in vertical equilibrium, but remained thin, 
a different Mach number relation is to be used (C90b): 

\2 



(M +7 + ^) 2 (M_ 7 + & 



(5a) 



2 + ( 7 -l)M 2 2 + ( 7 -l)M 2 ' 

In this case, the postshock Mach number for a strong shock is M + = 0.45, and the 
proton temperature in the postshock region is, 

T p ~ 0.96 x lO 12 ^" 1 K (56) 

For a spherical flow with 7 ~ 5/3, shocks can form if the accreted sub-Keplerian matter 
is already supersonic (this is true for models in vertical equilibrium if 7 ^ 1.5, C95a) 
but for a thinner flow even with 7 = 5/3 there are two saddle type sonic points and 
shocks could form. Subsequent behavior of the protons and electrons will depend on 
the heating and cooling processes. As the matter accretes, it is compressed due to 
geometrical effects and becomes hotter. At the same time, the electrons lose energy 
due to bremsstrahlung and Comptonization of the soft photons from the disk. They also 
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gain energy from the protons due to the Coulomb interactions, and cool the protons in 
the process. The energy equation which the protons and electrons obey in the postshock 
region is: 



where e is given by, 



|-(e+-) + (r-A)=0, (6) 
or p 



- - 2 P - 1 Jl 
e ~2 Vr + p( 1 -l) 2(r-l) + 2r 2 ' {) 



and T and A are the heating and cooling terms respectively. 

In the equation which governs the proton energy, we consider 7 = 5/3 and no 
heating. We consider the coolings to be due to inverse bremsstrahlung, and the transfer 
of energy to the electrons by Coulomb coupling. In the equation which governs the elec- 
tron energy, we use heating due to transfer of energy from protons through Coulomb 
coupling and cooling due to bremsstrahlung and Comptonization. If the electron tem- 
perature T e is high enough, kT e > m e c 2 (where, m e is the mass of an electron), we 
use 7 = 4/3, otherwise we use 7 = 5/3. Synchrotron losses could be trivially added to 
our equations. Since this requires a knowledge of unknown magnetic field contents, we 
ignore its effects for the time being. The temperature at the shock is not high enough to 
induce a significant pair production, so we ignore these effects. Within the convergent 
flow regime r < 2 — 3r g , in optically thick flow, pair production will take place which 
may be responsible for the annihilation lines observed in soft state. 

2.3 Radiative Processes 

The bremsstrahlung and Coulomb cooling rates are taken from the standard texts 
(e.g., Lang 1980). Correct Comptonization rates have been computed very recently 
(T94; TL95) and we discuss them here for the sake of completeness. We provide new 
results on the computation of reflection of the hard radiation from the cooler disk. 

2.3.1 Comptonization 
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Comptonization is the problem of energy exchange in the scattering of photons off 
the electrons. In the postshock plasma, as the accretion rate increases and the Thomson 
opacity, 

rr s 

Tt = cr T n e dr, (8) 

approaches unity, the escaping photons gain energy due to repeated scattering with the 
hot electrons. For a thermal non-relativistic electron distribution with a temperature 
T e , the average energy exchange per scattering is given by (hi/, kT e <C m e c 2 ), 



<Au> 4kT e - hu 



v m e c 2 



(9) 



When hu <C kT e , photons gain energy due to the Doppler effect and when hu 3> kT e , 
photons lose energy because of the recoil effect. As the radiation passes through the 
medium, the probability of repeated scattering by the same photon decreases exponen- 
tially although the corresponding gain in energy is exponentially higher. A balance of 
these two factors yields a power law distribution of the energy density: 

F v oc u~ a . (10) 

In the limit of high energies (hu 3> kT e ) the exponential hard tail is formed as a result 
of the recoil effect (the photons are unable to gain more energy than electrons actually 
have: hu oc kT e ). 

It is shown in TL95 that the power laws are the exact solutions of the radia- 
tive transfer kinetic equations taking into account the Doppler effect only. The related 
spectral indices are derived for an wide ranges of the optical depths and electron tem- 
peratures of the plasma cloud from the following transcendental equation as follows: 

a ~ ln[l + (a + 3)0/(1 + 0) + 4dl /a e 2 } ' (U) 

where 

3[(a + 3)a + 4]r(2a + 2) 

an (a) = ; — — ; , 

V ; (a + 3)(« + 2) 2 

and = kT e /m e c 2 which represents the appropriately weighted average of the electron 
temperature in the postshock region is computed from the model. Here, (5 is obtained 
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from the eigenvalue C = exp(— (3) of the corresponding eigenfunction for the appropriate 
radiative transfer problem. If we approximate the postshock region as a slab, for optical 
depths t > 0.1, 

7T , _ _i n„ s „.__q 7^-., , 1 

3r ' 



12(r + 2/3) 
and for very low optical depth r < 0.1, 

13 = hi 



-(1 - e - 1 - 35T0 ) + 0.45e- 3 - 7ro In 



(12a) 



1 



(126) 



T ln(1.53/T )_ 

If on the other hand, the postshock region is approximated as a spherical bulge, then 
for optical depths (r > 0.1), 



3(to + 2/3) 2 
and for very low optical depth (r < 0.11) 



(1 - e- a7ro ) + e~ 1AT0 



4, 

3r 



(13a) 



(3 = \n 



3r n 



(136) 



In Hua & Titarchuk (1995, hereafter HT95) the Comptonization spectra are given for 
the optically thin (r < 1) and thick (t > 1) regimes. There is a good approximate 
formula combining both regimes (HT95 Eq. [9]) which we use in our spectral calcula- 
tions. 



From the numerical simulation of the three dimensional flows (MLC94), it is ev- 
ident that the postshock region behaves more like a thick accretion disk (Rees et al. 
1982; Paczyhski & Wiita 1980). Thus, it is neither a slab nor a sphere, but approxi- 
mately a torus. Therefore, one expects the spectral index to be somewhat in between 
these two limits. However, as one goes inward, the optical depth rises and the contri- 
bution to the external hard radiation field diminishes. Therefore, either of these two 
assumptions should give sufficiently accurate results. The results we provide in the 
next section assumes the geometry of the hot region to be spherical. This assumption 
is made only to obtain an appropriately averaged electron temperature to compute the 
emergent hard spectra. 
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For our self-consistent calculation of the temperature distribution of the shock 
region we use the expression of the Comptonization enhancement factor E Comp as a 
function of the spectral index a and the effective low- frequency (temperature T r ) di- 
mensionless photon 

q XQ {u)xo~ l , if a < 1; 

if a > 1 . 



Ecomp ' 



«(« + 3)(l-ft£xr 1 ) ;f _ (14) 



(a + 4) (a - 1) 



where, 



_ q(q + 3)r(a + 4)r(a)r(l - a) , 
?xo W- r(2a + 4) [i X ° 1 

This formula, obtained by integration of the emergent spectrum over the photon energy 
covers asymptotic forms for both the regimes (a < 1 and a > 1) which are presented 
in ST80, Sunyaev & Titarchuk (1985, hereafter ST85). 



2.3.2 Computation of the Reflected Component 



Because a great deal of importance has been given in recent literature on the 
possibility of observation of the reflected component (e.g., Zdziarski et al., 1990; Done 
et al. 1992; Haardt et al., 1993), we address this issue here in detail. We find angle 
and energy dependent analytical solution (for E < m e c 2 ) which is applicable when 
the reflector (disk) is cold ( <, 10 5 K). For a hotter disk, the reflection is more efficient 
because the photoelectric absorption process is considerably weak (Kallman & Krolik, 
1986). We have incorporated this effect by neglecting the photoelectric absorption term 
(eq. 19-21 below) for T £ 10 5 K. Our computation of the equivalent width in terms of 
a cold disk reflection only gives an upper limit, since it is smaller for any hotter disk. 

In computation of the reflected hard component, we note that the reflection spec- 
trum from the disk is a result of scattering of the hard radiations from the shock by the 
relatively cold matter of the disk. They suffer photoelectric absorption at energies less 
than 10 keV with emission of the strong K a line and down-scattering (recoil effect) at 
energies higher than 50 keV. Thus, the reflection spectrum has a maximum at an energy 
around 15 — 20 keV. In order to estimate the reflected flux one requires a knowledge of 
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the fraction f sd of the shock hard radiation intercepted by the disk and the reflection 
properties of the disk. In the very nonrelativistic energy range < 30keV, photons scatter 
off electrons almost coherently and consequently the results of the classical reflection 
problem (see e.g., Chandrasekhar 1960) are applicable, i.e., the intensity reflected from 
the disk at angle ^ = cos _1 fi (with respect to the disk normal) when the incoming hard 
flux ttFh falls onto the disk at the angle 6> = cos _1 yU is as follows 

7 ^°) = — ^ • (15) 

Here, 

\ u = {[1 + [Y + Y ■ Q(E - 7.11keV)](7.8 keV/^) 3 }" 1 , (16) 

is the photon scattering probability in the presence of photoelectric absorption, H(fi, \ v ) 
is the well known H-function which changes within a range of 1—3 or even less depending 
on the photon scattering probability (see e.g., Chandrasekhar, 1960), Y and Y are the 
abundances of elements (in units of the cosmic abundance) with a charge Z < 26 and 
the iron abundance, respectively. Q(x) is the Theta- function (0 for x < and 1 for 
x > 0). At the hard tail, the intensity can be presented as a product of the three 
factors, namely /(/x, /x ) — F H f sd A u (ii, /io), where A u (/i, /i ) is the albedo of the cooler 
disk. 

Computation of Albedo 

In the computation of the albedo A u (r, /j, /x ) of the disk, same two processes 
are to be taken into account: the recoil effect and the photoelectric absorption. In 
the presence of the recoil effect alone, the incident radiation of frequency v will be 
completely absorbed, if (ST80, Titarchuk 1987, hereafter T87) 

Au hv 



r 2 = zrj ~ 1 (17) 



2 

v m e c z 
i.e., 

to - (is) 
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The albedo A u would be 1 — 1/t = 1 — \fz. In the presence of photoelectric absorption, 
similarly Tq(1 — A^) ~1 and A u = l — l/r = l — where, 

[y + y -e(E-7.iikeV)](^/^) 3 
" l + ir + ro-e^-T.iikeV)]^/^) 3 

(z* = 7.8 keV/511keV) is the probability of the photoelectric absorption (see, Eq.16). 

In the presence of both recoil and photoelectric absorption, and taking into ac- 
count the angular dependence of the incoming and outgoing radiations, we (see T87) 
obtain a more accurate expression for albedo as a function of the frequency of the 
incoming radiation and the angle of incidence 9 as, 

A v = 1 - 0(/i O )A, (19) 

where 

0M = — 1 -r== t 1 - ^r(* + A ')( ln ^ + 133 - i-458/ir 2 )], (20) 

1 + ^3(1 - K)fi 4 

A = (1 - \ u ) 1/2 = [v^exp(^)erfc (^) + 5 l J\ (21) 

z \/Z 



and fi = cosip where ip is the angle which the incident radiation subtends with the local 
normal to the disk. The second term is negligible for T £ 10 5 K. 



2.3.2.2 The analytical solution of the reflection component 



The problem of the reflection from the cold material is related to the radiative 
transfer kinetic equation which takes into account the photoelectric absorption and 
the recoil effect (see e.g. Basko, Sunyaev & Titarchuk 1974, Illarionov et al. 1979, 
T87, Grebenev & Sunyaev 1987, Magdziarz & Zdziarski 1995). We take into account 
contribution of the multiple scattering, using Fokker-Planck (diffusion) approximation 
and the first scattering exactly. This is a standard approach used in the classical 
radiative transfer theory for calculation of the coherent multiple scattering (e.g. Sobolev 
1975). 
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The illumination of the semi-infinite atmosphere (disk) by the parallel flux irF(z) 
with inclination angle with respect to the disk normal 9 = cos _1 // produce the expo- 
nential primary photon distribution throughout the atmosphere as follows 

f(T,z) = ^^eM-r/Vo), (22) 

where r is the Thomson optical depth calculated from the upper disk surface. The 
first scattering component of the emergent spectral intensity observed at the angle 
6 = cos -1 /! with respect to the disk normal is calculated by the integration of the 
primary photon distribution with the escape probability exp(— r//x)//x along the line of 
sight (see, e.g., Chandrasekhar 1960) 

h{H,Ho,z) = — ■ (23) 

4 fj, + Ho 

The multiple scattering component is produced by the integration of the Fokker-Planck 
(Kompaneets) equation and by integration of the product of its solution and the escape 
probability exp(— r///)//i along the line of sight. 

The Fokker-Planck (Kompaneets) equation reads (T87, GS87) 

1 d / / \ o \ 1 d 2 n a A \ v F(z) , . . 

"> + 3^)9^ - V T n = -^eM-rM- (24) 

where n(z, r) is the photon occupation number, 

— = 1 - X u ^ [Y + Y ■ Q(E - 7.11keV)](7.8 keV/hu) 3 , 

(TT 

<j)(z) = (1 + 2.8^ - 0.442 2 )- 1 , (25) 

'W = 1 + 4 .fc+l.lz» - (26) 
The solution of this equation with the appropriate boundary condition at To = 
(dn/dr - 3/2 = 0) is given by the integral (T87, GS87) 

J ^ z ^ o) = MzjL exp [-L w^ dt ) 



Zl 



-d Zl / °exp(-e//io)P[r, M (z,z 1 ),e]^. (27) 
Jo 
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where J(r, z, /x ) = z 3 n(r, z, /i ) is the average intensity, 

a/3 

P(r, u, = -Wexp[-3(r - 2 /4«] + exp[-3(r + if/Au] 

ZWTTU 



-3/ exp[-3(r + e + ^) 2 /4M-3r7/2]cir/}, (28) 

JO 

u(z, Zl ) = -- — + 7.4 ln(zi/z) + 13.54(z - (29) 

2 Z\ 



[ Z1 -^M-d* = 3.375 • l(T 6 [y + F • ©(# - 7.11keV)l 
J 2 rj{t)a T 

•[0.25(^ 4 - ^r 4 ) + 1.53(z -3 - ^r 3 ) + 0.55(^ 2 - z7/ 2 )]. (30) 

Evaluation of the internal integral in equation (27) over £ followed by the integration 
(over r) of the product of the average intensity J(r, z, /i ) and the escape probability 
exp(— t /[/,)/[/, along the line of sight gives us the multiple scattering component of the 
emergent spectral intensity 

Im ^^ z) = Mz)L exp [-l ^ dt ) 

^—^- G\p, no, u(z, zi)/3]rfzi. (31) 

Z\ 

where 

G ^^ t] = (^ - jflfo - 2/3) (/i - 2/3) ^ ' 4/9V(t ' 

-/x(/i 2 - 4/9)^(t, //) + 2/3(^ 2 - /i 2 )v9(t, 2/3)], (32) 

where t = u/3 and 

<p(t, ji) = erfc(\/t//z) exp(t//i 2 ). (33) 

Thus, the final expression of the emergent spectral intensity is a sum of two components 
(eqs. [23] and [31]) 

I(H, Ho, z) = h(n, n , z) + I m (n, no, z). (34) 



2.3.3 Spectral Properties of the Converging Inflow 



As the optical depth of the flow becomes larger than 1, the momentum deposition 
by the bulk motion on the cooler soft photons becomes important (e.g. Blandford & 
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Payne 1981; Lyubarskii & Sunyaev, 1982). In the soft state, when the postshock region 
becomes cool and optically thick, the spectral properties would be dictated by this 
converging inflow. 

The parameter which shows the importance of the bulk Comptonization effects 
relative to thermal Comptonization is 5 = l/9m. For kT e <C m e c 2 and m > 1, 5 >> 1 
and this leads to an asymptotic relation for the spectral index a 

a = 2A 2 - 3. 

Therefore, as long as 5 1, bulk Comptonization effects dominate and the values of a 
are indepedent of the temperature of the electrons. For details of derivations of these 
parameters, see TMK96. 

2.3.3.1 Spectral Index of Comptonization by Converging Inflow 

In order to understand the origin of the hard component in converging inflow, 
one requires to solve an equation of the spectral energy flux F(r, v) with the boundary 
conditions of F — > at r — > oo, and F = — 0.5nx 3 at the inner boundary: r — > 1 (half 
the flux is absorbed by the hole). This boundary condition along with the radiative 
transfer equation imply that eigenvalue A 2 is the root of the equation 

~ ^^t 6 ) $(-A 2 + 5/2, 7/2, r b ) + h -^ n ^{-\ 2 + 7/2, 9/2, r b ) = 0, 

where = 1.5m and $ is the well known Kummer's function (e.g. Abramowicz & 
Stegun, 1964). In the limit of the large m, the first root of this equation A 2 = A 2 2.25 
(TMK96) and thus a — 1.5. Though this slope is obtained from Newtonian consider- 
ations, similar result is likely to remain valid even in general relativistic computation. 
This will be discussed elsewhere. 

2.3.3.2 Luminosity and Effectiveness of Comptonization through Converging inflow 
The relation between the luminosity of the low frequency photon source L = 
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J °° 5(x — x )dx = 1 and that of the source of hard radiation photons subjected to the 
converging inflow Comptonization (x = hu/kT e ), 



/ I{x)dx (35) 
Jo 



is obtained using the integration technique developed in ST80, LS82, ST85, 



where £ = a + 4 + 5. 



1 '1-- — - ir 1 ], fora>l 



T{t)T{a)T{\-a){\-x\-«)xr X 

r(a + o 



(36) 



for a < 1 



This is a generalization of ST80 results for the converging inflow case. Eq. (36) 
also has a continuous transition through spectral index a = 1 (see eq. 14). 

The important conclusion from the above exercise is that the low-frequency source 
flux is amplified only by a factor of 3 (a ~ 1.5) due to Comptonization of the converging 
inflow into black hole. Thus the component is very weak and mainly visible in the soft 
state when either the postshock region is cold or when the shock is absent. 



2.4 Solution Procedure 



A fully self-consistent solution of the system shown in Fig. 1 involves the following 
considerations: 



(1) The postshock region is hotter and thicker: h(r s ) ~ a + (r s )r^ 2 , a + being the sound 
speed in the postshock region. It allows this region to intercept a fraction f^ s of the 
soft photons from the standard disk component. We compute self-consistently this 
intercepted soft component from each radial distance of the disk and integrate over 
the disk to obtain fa s . This intercepted radiation is Comptonized in the postshock 
disk. From the two-temperature equations (eq. 6) we obtain electron temperature as a 
function of the optical depth. We then obtain an average electron temperature and the 
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enhancement factor following the prescription of ST85 and Titarchuk (1988, hereafter 
T88). In reality, we approximate the numerical averaging curves of ST85 and T88 by 
the following accurate analytical function (valid for spherical geometry): 

9(r) = (1 - V(-+ 2 ))co^(l -L) + | e -(*.+2) 

and obtain the average electron temperature from, 

T( ^T e {r)g\r){r -rfdr 

Here, r is the total optical depth of the postshock region. 

(2) After the computation of the Comptonization process, we calculate the fraction f sd 
of this hard radiation that is intercepted by the standard disk. The interception is 
computed at each radial distance of the disk and then integrated over the whole disk. 
The rest (1 — f s d) of the hard radiation is directly radiated away to observers at infinity. 

(3) We compute the albedo A u (r) of the preshock standard disk which determines the 
fraction of the intercepted flux scattered away at each radius of the disk. The rest 
B„(r) = 1 — A v (r) is assumed to be absorbed by the disk and is reradiated. 

(4) In the zeroth order calculation, we assume that the equatorial Keplerian disk emits 
exactly as a standard Shakura-Sunyaev disk for r > r s . After it absorbs a fraction 
fsd&v{ r ) of the hard-radiation and reprocesses it, the temperature of the disk is com- 
puted at each radial distance. 

We now give a simple argument why our algorithm (1-4) should converge: Let Lj = 
fdsLss denote the fraction of the Shakura-Sunyaev (1973) disk luminosity L ss of the soft 
radiation intercepted by the bulge of the shock, Ec om p — Lh / Lj denote the enhancement 
factor of this radiation due to cooling of the electrons through inverse Comptonization, 
and f sd (~ 0.25 for a spherical bulge) denote the fraction of L H intercepted back by 
the disk. The soft component observed from a disk around the black hole candidate 
is given by equation (2). The enhancement factor E Co mp ~ 3(T e /3T d ) 1 ~ a ~ 10 — 30 
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(for a ~ 0.7 — 0.8) because, typically, the electron temperature T e ~ 50 keV and the 
disk temperature T d ~ 5 eV for parameters of active galaxies and T e ~ 150 keV and 
T d ~ 100 eV for stellar black hole candidates. Hence, we easily achieve a convergence, 
fdsEcompfsdBu < 1 since f sd ~ 0.25, f ds ~ 0.05, and B v ~ 0.5. 

In each run, we repeat steps (1) through (4) until the temperature distribution 
and the spectral index a converges. Typically, 5 to 10 iterations are enough to achieve 
convergence. 

3. RESULTS AND INTERPRETATIONS 

In each case we compute, we choose a set of four quantities, namely, the mass of 
the central black hole M, the shock location r s (instead of angular momentum of the 
halo), the accretion rate of the disk m d = M d /M Edd and the accretion rate of the halo 
rrih = Mh/MEdd- We shall present a series of runs for a galactic black hole candidate, 
M = 5M , and a series of runs when a black hole is massive, M = 1O 7 M . In both the 
cases, we shall consider parameters suitable for a non- rotating black hole (r s = 10r s , 
Ti = 3r g ). By choosing two independent accretion rates, we thus do away with the 
unknown viscosity parameter. 

In Fig. 2, we show the temperature distribution inside the disk and the postshock 
region when m d = 0.001 (solid line), 0.01 (long-dashed line), 0.1 (short-dashed line) and 
1 (dotted line), fn h = 1 is chosen throughout. We choose M = 5M . As the accretion 
rate of the Keplerian disk is increased, its temperature is increased. The postshock 
region, on the other hand, is cooled, since the number of soft photons intercepted by 
the postshock region is increased. The more the electrons cool through the inverse 
Compton process, the more the protons supply them energy through Coulomb coupling. 
Eventually, for m d = 1, the runaway process is seen where both the electrons and protons 
become cold catastrophically. Subsequently, the cooler but denser flow cools through the 
bremsstrahlung process alone. Exactly similar behavior is seen when the computation 
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around a massive black hole is also carried out. 



If Comptonization is the dominant mechanism for cooling and the densities and 

1 3 

velocities vary as in a spherical flow, i.e., v oc r~2 and p oc r^2, one can obtain an 
analytical expression for the electron temperature from the following equation (see, eq. 
[6]): 

1 <U " + --C Com y/>~0, (37) 



T e dr r 

where Ccomp (obtainable easily from eq. [6]) is a monotonically increasing function of 
optical depth and the spectral index is assumed to be constant for a given case. The 
solution of the above equation is 

Te = ^ e <Wr 3/2 -. 3/2 ). (38) 

This shows that if Ccomp > r~ 3 ^ 2 , the cooling due to Comptonization overcomes com- 
pressional (geometric) heating and T e drops as the flow approaches the black hole, and 
the convergent flow regime begins. Our numerical integration of equation (6) requires 
a somewhat higher amount of soft photons for cooling because of other heating effects, 
such as energy gain from the protons due to Coulomb process. The catastrophic cooling 
of the postshock region in soft states in important, since it introduces seed soft photons 
within the 'convergernt flow' regime and produce the weaker hard component of slope 
~ 1.5. 



In Figure 3a, we present the variation of the energy spectral index a (observed 
slope in the 2 — 50keV region) when a black hole of mass M = 5M Q is chosen. The 
abscissa is the logarithmic mass accretion rate of the Keplerian disk component (m^) 
and the different curves are marked by the mass accretion rates of the halo component 
(rhh) ■ We compared the results from the postshock flow as well as the convergent inflow. 
The dashed curves are drawn where the spectra could be contaminated by both effects. 
In Fig. 3b, we show mean electron temperature. The solid curve is for rhh = 1.0 and the 
long-dashed and short-dashed curves are for m h = 2.0 and 0.5 respectively. For a given 
halo rate, a transition between states can be achieved by a change in rh d . Qualitatively, 
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the cooling rate due to Comptonization for unit optical depth is proportional to the 
electron temperature and the local energy density, i.e., L H /rh oc T e (L H + Lj). Hence, 
ThT e oc (1 + LijLuY 1 ". In the hard state, Lj <C L H and r h T e ~ constant. When 
Th £ 1, this leads to a ~ constant (T94, TL95, see also, Pietrini & Krolik, 1995). This 
is what we observe as well: a is insensitive to the rhd when m d 0.1. In the soft state, 
Lj » L H and T e ~ t^t - - Smaller T e leads to a very high a, as observed. For a 
smaller rhh, the optical depth becomes small and the temperature high, causing the 
index to become especially sensitive to both and rhh- When accretion rate of the 
disk is very high, the spectral index asymptotically becomes 1.5, as is observed in many 
systems (e.g. Sunyaev et al. 1994, Miyamoto et al. 1991, Parmar et al. 1993, Wilson 
& Rothschild 1983). Successful fits of soft spectra of black hole candidates GS1124- 
68 (Ebisawa et al., submitted) and GRS 1009-45 (Titarchuk et al., submitted) have 
been obtained using our model. These fits are therefore consistent with a freely-falling, 
sub-Keplerian component close to a black hole. 

In Figures 4(a-b), we present the same quantities computed for a massive black 
hole (M = 1O 7 M ) candidate. The results are similar, though the average temperature 
is smaller and the index a is more sensitive to the disk rate when m d ^ 1 and fn h ^ 1. 
Here, the rapid rise of a is due to the abrupt cooling (see Fig. 2) of the postshock 
region within a narrow layer. 

Figure 5 shows the contribution of the component spectra for a particular case 
with m d = 0.1, m h = 1.0, M = 5M Q . For illustration, the observation angle 9 is chosen 
so that [i = cosQ = 0.4. The long dashed curve represents the soft spectrum from the 
Keplerian disk, while the short dashed curve is the contribution to the soft spectrum 
from the reprocessed and intercepted hard radiation by the Keplerian disk (second 
term in eq. [2]). The dotted curve is the contribution to the hard radiation from the 
postshock region. The dash-dotted curve denotes the reflected hard radiation from the 
Keplerian disk. The sum of these contributions is depicted by the solid curve. 
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Figure 6 shows a comparison of the four runs for the spectra around a black hole 
of mass 5M and m h — 1.0. The disk accretion rates (fn d ) are 0.001 (solid line), 0.01 
(long-dashed line), 0.1 (short-dashed line) and 1.0 (dotted line) respectively. As shown 
in Figs. 3a, 3b, 4a, 4b, with the increase of the disk accretion rate, the temperature 
T e of the electrons is reduced and the energy index a is increased. The luminosity 
and the peak frequency of the soft component go up monotonically with m^. The 
hard component shows a 'pivoting' property: the intensity of the hard component rises 
with m d below u ~ iq 18-19 Hz (corresponding energy ~ 5 — 50 keV) but it falls at 
higher energies since the electrons become cooler. In the soft state, the postshock 
region becomes cooler and produces a weak hard component of slope (a ~ 1.5) due 
to the convergent flow as shown by the dash-dotted curve. In Fig. 7 we show this 
'power-law' feature with the components separated. The long-dashed curve is the soft- 
luminosity from the Keplerian accretion disk and the short-dashed curve is due to the 
Comptonization of the intercepted (~ 5%) soft photons in the convergent flow. This 
power-law component is possibly the only true signature of black hole and is observed 
in black hole candidate spectra (e.g. Sunyaev et al. 1994). In the neutron star case, 
the flow need not pass through the inner sonic point (C89, C90a,b) if it is not compact 
enough and the flow is subsonic at the surface. Thus, no significant bulk momentum is 
added to the soft photons to create this hard component. Indeed, in this case, turbulent 
mixing of halo and disk components in the postshock region becomes important, which 
results in a hard component of slope < 1 and electron temperature T e ~ 12 keV which 
is observed (Claret et al. 1994). 

If the angular momentum of the overall disk, or even the sub-Keplerian compo- 
nent, is quite low, / <C l ms , the disk need not have a shock at all. Shock free advective 
disk (orginating out of Keplerian component, see, Appendix) can also intercept soft 
photons from the Keplerian disk and show similar behaviors mentioned here. It is un- 
clear whether X-ray emitted in this process would be efficient. As the viscosity of the 
disk varies, the sub-Keplerian flow would meet the Keplerian disk at varying distance 
(Appendix, C95b) and could mimic the hard-to-soft state transition. However, in this 
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case soft and hard components would always be correlated which is not observed. 

Increasing the halo rate m h increases the optical depth of the postshock region 
thereby increasing the hard component luminosity The spectrum is hardened due 
to the saturation of Comptonization (see, Figs. 3a and 3b). Figures 8a and 8b show 
the comparison of electron temperatures (in K) and the spectra, respectively, as fn h is 
varied while ma — 0.1 is kept fixed throughout. Solid, long-dashed, and short-dashed 
curves are drawn for rhh = 1,2 and 0.5 respectively. Contrary to what is observed in 
Fig. 6, the hardening of the spectra in the present circumstance is accompanied by no 
significant change in the soft luminosity and an actual increase in the hard luminosity. 

The above computations use the shock location to be r s ~ 10r g and the inner 
edge at 3r g . These are the typical length scales for accretion of matter with marginally 
stable angular momentum (C89, C90a,b). An increase in angular momentum of the halo 
increases the location of the shock, cools the average postshock electron temperature 
and raises the average spectral index. It also suppresses the soft luminosity relative to 
the hard luminosity. Similarly, a study of the behavior in the Kerr geometry with, say, 
r s = 5r g and r\ = 1.5r g with similar disk and halo parameters indicates a general increase 
of the average halo temperature and an average decrease of the spectral index. 

The above mentioned properties from our model are summarized (see also, CTKE95) 
in Table 1 where the correlation (arrows pointing up) or anticorrelation (arrows pointing 
down) of the observable quantities (luminosities L x ^ in X-ray and 7-ray regions) with 
the input accretion rates are shown. Smaller arrows represent a weaker correlation. 
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TABLE 1 

Variations of the Spectral Properties with Accretion Rates 





Input 






Output 

i 


Row 

& 

Column 


(1) 


(2) 


L s 


(1) 


T 


T 


L x ,a 


(2) 


TT a 




L 7 , a 


(3) 


IT* 





a dependence is weaker for rhd £ 0.1 

b dependence is weaker for rhh ~ 1 

c «i ~ 1 — > 1.5, L x /L s < 10" 3 in the convergent flow regime (see, Fig. 3) 

4. COMPARISON OF THE PREDICTED BEHAVIOR WITH 
OBSERVATIONS 

So far, we have delineated the spectral properties of the accretion disks which 
are formed out of Keplerian and sub-Keplerian matter. We can now compare these 
predictions with the observed facts. Galactic black hole candidates, such as GX339-4, 
GS2023+338 and GS1 124-68 show almost constant spectral slopes over two decades of 
the luminosity variation. (TAN; Ueda et al. 1994, hereafter U94; E94). The constancy 
of the slope is explained by the lower left corner of Fig. 3a. The spectral evolution 
of the black hole candidate novae, GS2000+25, GS2023+338, and GS1124-68 are very 
different (TAN; E94). These differences could be accounted for by variations of the 
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properties of the halo (Figs. 8a and 8b). For instance, GS2023+338 is always in the 
hard state throughout the outburst, suggesting a high m h > 1 but lower m d < 0.01—0.1. 
The suppression of the soft component could also be due to a distant shock. GS2000+25, 
on the other hand, remained in the soft state during the outburst, suggesting a high 
m d ~ 1.0, but a low m h £ 0.1 — 0.5. In the rising phase of GS1124-68 (E94), the 
increase of L H below ~ 10 keV is accompanied by its decrease above ~ 10 keV. This is 
the 'pivoting' property in the hard spectra we referred to (Fig. 6). This object, along 
with other black hole candidates such as Cyg X-l and GX339-4 are observed in two 
distinctly different states. Generally speaking, in the hard state, a ~ 0.7 remains almost 
constant but in the soft state the hard component is very weak with index a ~ 1.0 — 1.5, 
again remaining constant during a the period in which L s varied by as much as two 
orders of magnitude (TAN, E93, E94). In the soft state, m d is very high and the 
postshock optically thick region produces a spectral index of ~ 1.5 as it would in a 
convergent flow (Fig. 6). In GS2023+338, the tendency of the anticorrelation between 
L x and L 7 is also seen (TAN, Sunyaev et al., 1991), but they do not actually cross, 
possibly because both m d and rhh change together. The hard component in galactic 
black hole candidates is generally found to be more variable than the soft component 
and the variations are found to be on time scales of days to months (TAN, 192, E94). 
This is possibly due to independent variations in accretion rates as we have chosen in 
our model. 

Among the massive black hole candidates associated with AGNs, NGC 4-151 has 
been studied extensively (Yaqoob 1992; Yaqoob et al. 1993). NGC 4151 shows a to 
be generally stable though it is seen to rise with the hard luminosity in the 2 — 20keV 
range. This is consistent with our result of variation of m d (Fig. 6) where we show 
that below the pivoting, the hard luminosity rises with a. The a-Ln plot also showed 
some scatter. This can be understood from our Fig. 8b, where we show that when 
the halo rate is increased, the hard flux increases with the decrease of a. Because of 
longer time-scales, Seyferts do not show a change of states. We believe that possibly 
all Seyfert-ls belong to different regions of the same diagram (Figures 4a and 4b). A 
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spectral index of ~ 1.5 is also reported in MKN841 (Arnaud et al., 1986) whose origin 
could be the same as that of the galactic black holes discussed above. 

5. CONCLUDING REMARKS 

In this paper, we have presented the spectral properties of the accretion flows 
which are primarily sub-Keplerian or an admixture of Keplerian and sub-Keplerian at 
the outer boundary. Such matter produces an optically thick standard disk on the 
equatorial plane and the sub-Keplerian optically thin halo surrounding it produces a 
shock close to the black hole or a neutron star. We showed that the reprocessing of 
the soft-photons from the standard disk by the hotter, postshock region of the disk is 
sufficient to produce the observed hard radiation with correct behavior of the luminosity 
and the spectral slope. We showed that the the spectral index a ~ 0.7—0.8 could remain 
generally stable even when the accretion rate of the disk {rhd) is changed by more than 
2 orders of magnitude, as long as rhd ~ 0.1. This is seen in objects such as, GX339-4, 
GS2023+338 and GS1124-68 and possibly in some class of Seyferts. When the disk 
accretion rate is close to the Eddington rate, the sensitivity of a to the accretion rate 
goes up and the object switches to the soft state as the disk accretion rate is increased. 
We also showed that during the change of the state, the hard luminosity remains the 
same in some energy band, a property (pivoting) that is observed in GS1124-68. Usually, 
this location of the pivoting depends upon the accretion rate. A search for this property 
in a wide band (2 — 200keV) should be made in other black hole candidates, such as 
Cyg X-l and GX339-4 to verify our model. A similar property is expected in X-ray 
bursters, and other neutron star candidates as well. We also showed that in the soft 
state, the consideration of the convergent flow becomes applicable. We believe that 
the spectral slope of 1.5 that is seen in the soft state is the signature of a very large 
converging flow close to the black hole. For a neutron star accretion, such a component 
should clearly be absent. Traces of convergent inflow should be observable at very high 
energies even in hard state, since the halo rate remains high. 
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It may be possible that for stellar black holes, the two components of the disk, 
namely, the sub-Keplerian and the Keplerian, are of different origin. The sub-Keplerian 
component may be from the winds, whereas the Keplerian component is primarily from 
the surface of the companion. If they were of same origin (i.e. from a completely 
Keplerian or a sub-Keplerian disk), the variation at the source would be reflected in the 
variation of the hard and the soft components at two different times. This is because 
the halo is almost freely falling, whereas the disk is moving on a viscous timescale. 
Roughly speaking, for a similar halo and disk thickness variation with radius, the ratio 
of the time scales would be tdisk/thaio ~ l/cx v , where a v is the Shakura-Sunyaev viscosity 
parameter (Shakura & Sunyaev 1973) inside the disk. However, even this correlation 
may be generally difficult to detect, since, whereas the variation of the disk rate changes 
both the soft and the hard luminosities, the variation of halo rate changes primarily the 
hard component. Also, the variation of the hard flux and the spectral index are inversely 
correlated with these rates (Figs. 6 and 8b). Another source of complexity could be 
the entrainment and evaporation of some disk matter which joins the halo and the loss 
of some matter as winds and outflows. Details of these processes can be quantitatively 
understood through numerical simulations and will be discussed elsewhere. 

We have noted above that there could be two sources of strong winds in the 
vicinity of the disk. Numerical simulations show that centrifugal barrier pushes sub- 
Keplerian, inviscid flows along the direction of the jets (Hawley et al. 1984, 1985; 
Eggum, Coroniti, & Katz, 1987; MLC94, Ryu et al. 1995). Another source of winds 
could be the evaporation of matter close to the shock. Judging from the resemblance of 
the iron resonance line shapes (Piro, Yamauchi, & Matsuoka, 1990; Rajeev et al. 1994; 
Mushotzky et al. 1995) with the P-Cygni profiles and the profiles of down-scattered 
emission lines we believe that these lines could be generated from the winds themselves. 
These profiles are characterized by the prominent emission feature at the red wing and 
by possibly weaker absorption feature at the blue wing. We are more convinced by the 
fact that we find the upper limit of the equivalent width of the lines from the disk to be 
only about a few tens of electron volts, as opposed to the observed equivalent width of 
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several hundreds of electron volts. A back-of-the-envelope calculation for the ionization 
parameter (Kallman & McCray 1982) could be done using typical parameters for the 
wind. This turns out to be, £ = L H \7-iokev/n p r 2 ~ 10 4 (^ 2l )2(^)o.i where, n p is the 
number density of protons and m w is the outflow rate of the line driven wind in units of 
the Eddington rate. Here, a tenth of velocity of light (subscript 0.1) and a 50% of loss in 
halo matter (subscript 2) are assumed. The column density -/V = 3.5(m/ l )i(^) Ch i(^-)~ 1 
g cm~ 2 may be sufficient to produce the observed iron line in the outflow. The subscript 
1 of rhh indicates that one Eddington rate is assumed. The emergent spectrum of the 
outflow is a result of transmission of the injected hard radiation. Fig. 9 shows the fate 
of this hard component (with a typical a ~ 0.7 for illustration) when it passes through a 
wind of r ~ 3. The solid curve shows the emitted component from the postshock region 
and the long-dashed curve is the component down-scattered (ST80) by the outflow (for 
energy change by down scattering, see, eq. 17). The short-dashed curve is produced due 
to the absorption by heavy elements at the outer layer of the outflow (pronounced iron 
edge is seen). We do not include the line emissions. If the so-called "reflection bump" 
originates in the wind, as we suggest, the net effect must be a combination of the solid 
and the short-dashed curve weighted appropriately by the covering factor. In computing 
the effect of reflection, it is well known that only a finite layer of an otherwise semi- 
infinite target participatesin scattering. Thus, the continuum spectra from these two 
approaches (namely, from reflection and transmission) are indistinguishable. However, 
the equivalent width of the emitted lines should be higher in the outflow, because of a 
larger (by an order of magnitude or so) covering factor. The equivalent width of the 
lines and the strength of our "transmission bump" should be correlated as they both 
depend on the covering factor. In the rotating winds and outflows, splitted emission 
lines should be generated. Such observations are reported in Tanaka (1994) and Tanaka 
et al. (1995). High optical depth (r > 3) of the disk components close to the hole washes 
out any disk emission features and thus these lines can only be emitted from the winds. 
Down scattering mimics gravitational red-shifts of lines. Thus the line profiles in these 
two processes are again indistinguishable, apart from the sharp cut-off in the red wing 
when a disk emission is considered. 



33 



Combining results presented in Fig. 1, Fig. 3 and Fig. 10, we are now in a position 
to construct sequences of a nova outburst and its quiescent state (e.g., Cannizzo, 1993). 
In the quiescent state, entire disk has low viscosity and the Keplerian disk receeds very 
far away (Fig. 10). As the viscosity in the equatorial plane rises, Keplerian disk rate is 
increased, and the inner edge of the Keplerian component comes closer to the black hole. 
This increases soft photons to be Comptonized, thereby generating an outburst. This 
situation is to be compared with an atomic reactor whose energy output is controlled 
by insertion of the fuel rods (here, the Keplerian disk component) by varying degree. 

In our model we have made some approximations, such as the equations appropri- 
ate for conical geometry (instead of actual torus) to integrate both the Euler equations 
as well as the radiative transfer equations. Also, we assumed a fixed polytropic index 
rather than computing it self-consistently form the heating and cooling as the flow 
accretes. We have ignored diffusion of hard photons to the preshock region and subse- 
quant pre-heating of the electrons, as is the characteristics of the photo-hydrodynamical 
shocks (Riffert, 1988; Becker 1988). However, we do not believe that our basic conclu- 
sions will be affected by these approximations. 

It is to be noted that our model differs significantly from the Comptonized soft 
photon model of two-temperature hot disks discussed almost two decades ago (Shapiro, 
Lightman, & Eardley 1974; Shakura & Sunyaev 1976). In these models, the soft and 
the hard radiations originate from the same disk and therefore, the variations of the 
hard and the soft components should be always correlated, this does not appear to 
be observed. There are some other models in the literature (Wandel & Liang 1991; 
Haardt & Maraschi 1991) which required extra-components and/or extra parameters. 
Our present model is somewhat more general than the single component model (which 
included shock waves) of Chakrabarti & Wiita (1992) where only the optical and UV 
regions were considered. 

If the shock waves do exist, as we assumed in our model, there could be another 
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manifestation which is observable as well. Recently Molteni, Sponholz, & Chakrabarti 
(1996) have shown that when the accretion rates are such that the cooling time scale 
roughly matches the infall time scale, the shock oscillates with a time period comparable 
to the cooling time. Could this be the origin of quasi periodic oscillations (QPOs) seen 
in black hole candidates (e.g., Dotani 1992) such as GX339-4 and GS1124-68? Possibly, 
yes. Particularly noteworthy is the fact that by dynamically moving back and forth, the 
shocks modulate the outgoing hard radiation by as much as 10 to 15 percent, whereas 
the soft radiation emitted from the preshock disk is only marginally modulated. This 
is because of the fractional variation of the emitting area is larger for the postshock 
region. These behaviors have been observed in dwarf novae, such as SS Cygni, VW 
Hyi, and U Gem (see, e.g., Mauche, Raymond, & Mattei, 1995). Unlike the case of a 
dwarf novae outburst, where cooling could be due to bremsstrahlung or line effects, in 
the black hole candidates, the time scale is governed by Comptonization processes. 

Though we have primarily discussed shock waves around a galactic or extragalac- 
tic black holes, the entire computation remains unchanged when a neutron star is con- 
sidered instead. The typical size of a neutron star (10 km) is under 2.5r g and therefore 
our computation zone (3 — I0r g ) is outside that of a neutron star. For example, the best 
fit of the hard spectra of the burst source MXB 1728-34 by a Comptonization model 
(ST80) requires the parameters kT e = 12 keV and the optical depth r > 3.8 (Claret et 
al. 1994) which are easily achievable by our model with a disk and a postshock bound- 
ary layer. Though it may be generally difficult to distinguish the spectral properties 
of a galactic Schwarzschild black hole from that of a neutron star, we believe that the 
hard component due to the convergent flow (Fig. 6) is possible only in a black hole 
accretion since the inner boundary conditions are completely different (C89, C90a,b). 

Works of SKC and LGT are partially supported by the National Research Council 
and NASA grant No. NCC5-52, respectively. The authors thank D. Kazanas, K. 
Ebisawa, T. Yaqoob and T. Kallman for useful discussions. They also thank Paul 
Wiita for carefully reading the manuscript and suggesting improvements. 
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APPENDIX 



The globally complete viscous, isothermal solutions around a black hole are in 
C90a,b where sub-Keplerian flows become Keplerian at large distance. In this Ap- 
pendix, we show this even in the presence of general heating and cooling. 

We choose a simple, vertically averaged, axisymmetric, polytropic disk. The 
equations governing the flow are (C95a,b): 

(a) The radial momentum equation: 



dv IdP , l 2 Kep -l 2 



dr p dr r 3 

(b) The continuity equation: 



^L(prhv) = 0, (A.lb) 
dr 



(c) The azimuthal momentum equation: 



dl(r) 1 d^aPr 3 hdQ^_^ (A 1 ) 

dr prh dr Vlxep dr 



(d) The entropy equation: 



ds 

T,vT-f- = Q + -Q- (Aid) 
dr 

These equations are the generalization of viscous, transonic, isothermal disks for which 
the global solutions are already in the literature (C90a,b). Here lx e p and flxep are 
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the Keplerian angular momentum and Keplerian angular velocity respectively, £ is the 
density p vertically integrated, P is the total pressure, h is the vertical thickness of the 
disk at radial distance r, v is the radial velocity, s is the entropy density of the flow, Q + 
and Q~ are the heat gained and lost by the flow. We compute h(r) assuming the disk 
is in hydrostatic balance equation in vertical direction. l(r) is the angular momentum 
distribution of the disk matter. Here, we have chosen geometric units, thus r is units 
of Schwarzschild radius r g = 2GM/c 2 , l(r) is in units of 2GM/c, and velocities are 
in units of the velocity of light c, M being the mass of the central black hole. To 
mimic the black hole geometry (i.e., to compute the Keplerian quantities above) we use 
Paczyhski-Wiita (1980) potential. In eq. (A.lc), a v £ 1 in the above equation is the 
viscosity parameter of Shakura & Sunyaev (1973), which is widely used to describe the 
viscous stress: w re j, = —aP. Eqn. A.lc could be easily integrated to obtain the disk 
angular momentum distribution. l(r) — l(r in ) = a v ra 2 /v (C90a,b). Without any loss of 
generality, we assume Q_ to be a fraction of Q + while integrating above equations. Also, 
in order that the angular momentum remains continuous across shock waves, we choose 
total pressure (thermal plus ram) in the viscosity prescription (CM95). One needs to 
supply five quantities (over-determined system, since we are not explicitly computing 
the heating and coolings to illustrate our point here): the angular momentum at the 
inner edge l(ri n ), the location of the inner sonic point r c , the viscosity parameter a v , 
the polytropic index 7 = 4/3 and the cooling law Q_/Q + . 

Fig. A.l shows examples of the ratio of the disk angular momentum distribution 
to the Keplerian distribution (for r < 3r g , we keep /i^ ep (r) = /xep(3r 3 ) for stability 
reasons) as a function of distance from the black hole. The ratios show deviation from 
Keplerian due to advection, pressure and viscous effects. Three diverse cases have 
been chosen to illustrate our points. In Case A (marked 'A'), we choose l(r in ) = 1.88, 
r c = 2.2, a v = 0.005, 7 = 4/3 (radiation pressure dominated) and Q_ = Q + . The 
flow deviated from Keplerian disk at 7.5r g and even becomes super-Keplerian (ratio 
> 1) close to the black hole. In Case B (marked 'B'), we choose r c = 2.3, l(r in ) = 1.7, 
7 = 5/3 (gas pressure dominated), a v = 0.02 and Q_ = Q + . Here the flow deviates 
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from Keplerian at 90r g and always remained sub-Keplerian. In both the cases above, 
there is no shock formation. In Case C (marked 'C'), the flow passes through a shock at 
13.9r 9 (angular momentum distribution remained continuous) and remains completely 
sub-Keplerian after deviating from the Keplerian disk at 480r 9 . The parameters chosen 
are l(r in ) = 1.6, r c = 2.87, 7 = 4/3, a v = 0.05 and <5_ = 0.5Q+- It is easily shown that 
a decrease of viscosity, keeping other parameters fixed, increases the distance where 
Keplerian disk begins. This property is crucial in understanding our model of the 
generalized disk (Fig. 1) and the nova outbursts. 



38 



REFERENCES 



Abramowitz, M., & Stegun, LA. 1964, Handbook of Mathematical Functions, U.S. 
Dept. of Commerce, Washington. D.C. 

Abramowicz, M., Czerny, B., Lasota, J. P., & Szuskiewicz, E., 1988, ApJ, 332, 646 
Arnaud, K. et al. 1985, MNRAS, 217, 105 Basko, M.M. 1978, ApJ, 223, 268 
Basko, M.M., Sunyaev, R.A., & Titarchuk, L.G. 1974 A&A, 31, 249 
Becker, P., 1988, ApJ, 327, 772 

Blandford, R.D. & Payne, D.C 1981, MNRAS, 194, 1033 

Cannizzo, J.K., 1993 in Accretion Disks in Compact Stellar Systems, ed. J. Craig 
Wheeler (Singapore: World Scientific), 6. Chakrabarti, S.K. 1989, ApJ, 347, 365 (C89) 
Chakrabarti, S.K. 1990a, MNRAS, 243, 610 (C90a) 

Chakrabarti, S.K. 1990b, Theory of Transonic Astrophysical Flows (Singapore: World 
Scientific, Singapore, 1990) (C90b) 

Chakrabarti, S.K. 1993, Numerical Simulations in Astrophysics, Eds. J. Franco et al. 
(Cambridge Univ. Press: Cambridge) (C93) 

Chakrabarti, S.K. 1994, in Proceedings of the 17th Texas symposium, (New York 
Academy of Sciences, New York) 

Chakrabarti, S.K. 1995a, in Accretion Processes on Black Holes, Physics Reports (in 
press) (C95a) 

Chakrabarti, S.K. 1995b, ApJ Letters, submitted, (C95b) 

Chakrabarti, S.K., & Molteni, D. 1995, MNRAS, 272, 80 (CM95) 

Chakrabarti, S.K., Titarchuk, L., Kazanas, D., & Ebisawa, K. 1995, A&A Suppl. Ser., 

Proceedings of 3rd Compton Symposium (CTKE95) 

Chakrabarti, S.K., & Wiita, P.J. 1992, ApJ, 387, L21 

Chandrasekhar, S. 1960, Radiative Transfer, (New York: Dover) 

Chen, X. M., & Taam, R. 1993, ApJ, 412, 254 

Chen, X. M., Abramowicz, M., Lasota, J. P., Narayan, R. & Yi, I. 1995, ApJ, 443, L61 
Clavel, J. et al. 1990, MNRAS, 246, 668 
Claret, A. et al. 1994, ApJ, 423, 436 



39 



Done, C. et al. 1992, ApJ, 395, 375 

Dotani, Y, 1992 in Frontiers in X-ray Astronomy, (Tokyo: Universal Academy Press ), 
152 

Ebisawa, K. et al. 1993, ApJ, 403, 684 (E93) 

Ebisawa, K. et al. 1994, PASJ, 46, 375 (E94) 

Eggum, G.E., Coroniti, F.V., k Katz, J.I. 1987, ApJ, 323, 634 

Grebenev, S.A., & Sunyaev, R.A. 1987, Sov. Astron. Lett. 13, 438 

Haardt, F. et al., 1993, 411, L95 

Haardt, F., & Maraschi, L., 1991, ApJ, 380, L51 

Hawley, J.F., Smarr, L.L., & Wilson, J.R. 1984, ApJ, 277, 296 

Hawley, J.F., Smarr, L.L., & Wilson, J.R. 1985, ApJS, 55, 211 

Hua, X., & Titarchuk, L.G. 1995, ApJ, 449, 188 (HT95) 

Illarionov, A.F., Kallman, T., McCray, R., & Ross, R.R. 1979, ApJ. 228, 279 

Inoue, H. 1992, in Frontiers of X-ray Astronomy, ed. Y. Tanaka & K. Koyama, (Tokyo: 

University Press), 291 (192) 

Kallman, T.R., & Krolik, J. H. 1986, ApJ., 308, 805 

Kallman, T. R., & McCray, H. 1982, ApJ Supp. Ser, 50, 263 

Liang, E.P, & Thompson, K.A. 1980, ApJ, 240, 271 

Lang, K.R., 1980, Astrophysical Formula, (Springer Verlag, New York) 

Lybarskii, Yu.E., & Sunyaev, R.A. 1982, Soviet Astr. Lett. 8, 330 

Magdziarz, P. & Zdziarski, A.A., 1995, MNRAS, 273, 837 

Malkan, M. 1982, ApJ, 254, 22 

Mauche, C.W., Raymond, J.C., & Mattei, J. A. 1995, ApJ (in press) 
Molteni, D., Sponholz, H., & Chakrabarti, S.K. 1995, ApJ (submitted) 
Molteni, D., Lanzafame, G., & Chakrabarti, S.K. 1994, ApJ, 425, 161 (MLC94) 
Muchotrzeb, B., & Paczyhski, B. 1982, Acta Astron. 32, 1 

Miyamoto, S. et al. 1991, ApJ, 383, 784 Narayan, R., & Yi, I. 1994, ApJ, 428, L13 
Novikov, I., & K.S. Thorne. 1973. in: Black Holes, eds. C. DeWitt and B. DeWitt 
(Gordon and Breach, New York) 

Paczyhski B., & Bisnovatyi-Kogan, G. 1981, Acta Astron. 31, 283 



40 



Paczynski B., & Wiita, P.J. 1980, A&A, 88, 23 

Parmar, A.N. et al. 1993, A&A, 279, 179 Perola, G.C. et al. 1986 ApJ, 306, 508 

Peterson, B.M. et al. 1991 ApJ, 368, 119 

Pietrini, P., & Krolik, J.H. 1995, ApJ, in press 

Piro, L., Yamauchi, M., & Matsuoka, M. 1990, ApJ, 360, L35 

Rajeev, M.R., Chitnis, V.R., Rao, A.R. & Singh, K.P 1994, ApJ., 424, 376 

Rees, M.J., Begelman, M.C., Blandford, R.D. & Phinney, E.S. 1982, Nat, 295, 17 

Riffert, H. 1988. ApJ, 327, 760 

Ryu, D., Brown, G., Ostriker, J., & Loeb, A. 1995, ApJ, (in press) 
Shakura, N.I., & Sunyaev, R.A. 1973, A&A, 24, 337 
Shakura, N.I., & Sunyaev, R.A. 1976, MNRAS, 175, 613 

Shapiro, S.L., & Teukolsky, S.A., Black Holes, White Dwarfs and Neutron Stars, 1984 
(John Wiley & Sons: New York) 

Sobolev, V.V., 1975, Light Scattering in Planetery Atmosphere (Pargamon Press: Ox- 
ford) 

Sunyaev, R.A. et al., 1991, Astron. Lett., 17, 123 

Sunyaev, R.A. et al., 1994, Astron. Lett., 20, 777 

Sunyaev, R.A. & Titarchuk, L.G. 1980, A&A, 86, 121 (ST80) 

Sunyaev, R.A. & Titarchuk, L.G. 1985, A&A, 143, 374 (ST85) 

Sun, W.H. & Malkan, M. 1989, ApJ, 346, 68 

Shapiro, S.L., Lightman, A.P, & Eardley, D.M., 1974, ApJ, 204, 187 

Tanaka, Y., 1989, in the Proceedings of the 23rd ESLAB symposium, ed. J. Hunt & B. 

Battrick vol. 1, p.3 (Paris:ESA) (TAN) 

Tanaka, Y. 1994, in New Horizons of X-ray Astronomy, eds. F. Makino & T. Ohashi 
(Tokyo: Universal Academy Press) 
Tanaka, Y. et al. 1995, Nat., 375, 659 

Titarchuk, L.G. 1987, Soviet Astrofizika, 26, 57 (Astrophysics 1988, 26, 97) (T87) 
Titarchuk, L.G. 1988, Soviet Astrofizika, 29, 634 (Astrophysics 1988, 26, 97) (T88) 
Titarchuk, L.G. 1994, ApJ, 434, 570 (T94) 
Titarchuk, L.G., & Lyubarskij, Yu 1995, ApJ, in press (TL95) 



41 



Titarchuk, L.G., Mastichiadis, A. & Kylafis, N., 1996, ApJ, submitted (TMK96) 
Ueda, Y., Ebisawa, K., & Done, C, 1994, PASJ, 46, 107 
Wandel, A., & Liang, E.P., 1991, ApJ, 380, 84 

Wilson, C.K. & Rothschild, R.F. 1983, ApJ, 274, 717 Yaqoob, T., 1992, MNRAS, 258, 
198 

Yaqoob, T. et al, 1993, MNRAS, 262, 435 
Zdziarski, A.A. et al. ApJ, 363, LI 1990 



42 



FIGURE CAPTIONS 



Fig. 1: Schematic diagram of the accretion processes around a black hole. An optically 
thick, Keplerian disk which produces the soft component is surrounded by an optically 
thin sub-Keplerian halo which terminates in a standing shock close to the black hole. 
The postshock flow Comptonizes soft photons from the Keplerian disk and radiates 
them as the hard component. Iron line features may originate in the rotating winds. 

Fig. 2: Proton and electron temperatures (T p and T e in °K) inside the disk and the 
postshock region when m& = 0.001 (solid line ), 0.01 (long-dashed line), 0.1 (short- 
dashed line) and 1 (dotted line). Other parameters are rhh = 1 and M = 5M Q . At 
higher accretion rates, the postshock region cools faster and the convergent flow regime 
begins. 

Fig. 3(a-b): Variation of the (a) energy spectral index a (observed slope in the 2— 50keV 
region) and (b) the mean electron temperature (in keV) T e as functions of the disk and 
halo accretion rates. Notice the stability of a in the hard state and the transition 
of states when m d ~ 0.1 — 1.0. M — 5M is chosen. Halo rates are marked on the 
curve. In (a), computed spectral index in soft states due to convergent flow is also 
provided. Dashed curves indicate regions where both components could contribute. 
In (b), short-dashed, solid and long-dashed curves are for halo rates 0.5, 1.0 and 2.0 
respectively. 

Fig. 4(a-b): Same as Figs. 3a and 3b except that M = 1O 7 M is chosen. In this case, 
the average electron temperature is generally lower, and the spectral index generally 
higher. 

Fig. 5: Contributions of various components to the net spectral shape (solid). Long 
dashed, short dashed, dotted and dash-dotted curves are the contributions from the 
Shakura-Sunyaev disk r > r s , the reprocessed hard radiation by the Shakura-Sunyaev 
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disk, reprocessed soft-radiation by the postshock disk r < r s and the hard radiation 
reflected from the Shakura-Sunyaev disk along the observer (/j — OA). The contribution 
of the preshock halo which is orders of magnitude lower is ignored. Parameters are 
171(1 = 0.1, rhh = l, and M = 5 M . 

Fig. 6: Variation of the spectral shape as the accretion rate of the disk is varied. 
m d = 0.001 (solid line), 0.01 (long-dashed line), 0.1 (short-dashed line), and 1 (dotted 
line). Though the soft luminosity increases with rhd, the hard luminosity may show 
'pivoting' or crossing as intermediate energies. Note the correlation and anticorrelation 
of the hard flux and the spectral index before and after the pivoting point. The dash- 
dotted curve represents the hard component from convergent inflow near the black hole 
and has the characteristics of the slope ~ 1.5 in soft state. 

Fig. 7: Spectral shape of a Keplerian disk with a very high disk accretion rate. The 
postshock region cools immediately due to Comptonization and the cooler photons are 
subsequently Comptonized by momentum deposition from the converging bulk motion 
close to the black hole. The spectral slope is a ~ 1.5, an universal feature of black hole 
candidates in the soft state. 

Fig. 8(a-b): Variation of the (a) temperature (in °K) and (b) spectral shape as the 
accretion rate of the halo is changed. rfih = 1 (solid), rhh = 2 (long dashed) and m/ l = 0.5 
(short dashed). The hard flux goes up but the spectrum hardens with rhh while the 
disk temperature and the soft luminosity remain virtually unchanged. 

Fig. 9: Fate of a power-law hard component in presence of a wind of r ~ 3 with 
a hundred percent convering factor. Emitted (solid line), down-scattered component 
(long-dashed line) and its absorption by the cooler outer region (short-dashed line) 
are shown. This "transmission bump" from the wind is indistinguishable from the 
"reflection bump" from the disk, except for the higher equivalent width in the former 
case. 
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Fig. 10: Ratio of the disk angular momentum to Keplerian angular momentum in three 
different cases. In case marked 'A', the flow becomes sub-Keplerian before becoming 
super-Keplerian close to the black hole. In case marked 'B', the disk becomes sub- 
Keplerian and remained so before entering the black hole. In case marked £ C, the disk 
becomes Keplerian, then passed through a standing shock at r = 13.9r g before entering 
the black hole. 
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